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Abstract. We discuss properties of kinetic solutions of scalar conservation laws in the 
variational approach developed by Panov[7][H] and also Brenier[5]. Our main result shows 
that such solutions can be considered as curves in a suitable Hilbert space with tangents 
that are unique minimizers of an interaction functional. 



1. Introduction 

We consider a scalar conservation law 

«t + V a •/(«) = 0, (i,x)eR+xR n , (1) 

with the flux / : R — > R n . A function u(t, x) is called an entropy weak solution if for any 
convex entropy/entropy flux pairs (r](u(t, x)), q(u(t, x)) of the flux function /, 

Vt + V x -q<0, D'(R™ +1 ). (2) 

The existence of unique entropy solutions for (pQ) with the initial data u(t = 0) = uo 6 
L°°(M n ), as well as their stability, was obtained in Kruzhkov[5]- The solution can be 
described using the kinetic formulation, as proved by Lions-Perthame-Tadmor[6j. In this 
approach, u(t, x) is a weak entropy solution iff the kinetic density function 

Y ^x,v)= | Q v<u \ t ^ (3) 

verifies the transport equation 

Y t + f v (v) -V X Y= - 8 v m, P'(R™+ 2 ), (4) 

for some measure m 6 .M+(R™ +2 ). To be precise, in [6], the kinetic denisty x(t,x,i>) = 
H(v) — Y(t,x,v), where H = 1, v > 0, H = 0, v < 0, was used, but the result can be 
expressed through Y(t, x, v) as well. 

Condition (JU) can be equivalently expressed via a variational form: for any regular test 
functions Y(x,v), nondecreasing in v, it holds 



(Y - Y)(Y t + /„ • V,y) dxdv > 0, V'(R). (5) 
Indeed, the equivalence holds because, 

J f Y(Y t + /„ • V X Y) dxdv = dt f \u(t, x)\dx = 0, 
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if, for example, uq G L 1 (M n ), or periodic. An interesting kinetic formulation for weak 
solutions were obtained in Panov[71[8], by allowing generic nondecreasing in v functions 
Y(t,x,v) in (J5]), rather than functions with the range in {—1,0, 1}. 

Panov[7J introduced a class of strong measure- valued solutions on ([1]) as the set of 
parametrized probability measures i^t^x ^ • / ^+(-^)? fo-^ which the level sets of their distri- 
bution functions Y(t,x,v) := u ttX ((—(x>,v]) are the graphs of weak entropy solutions of 
©, i.e. 

u(t,x,X) = sup{v : Y(t,x,v) < A}, (6) 

is an entropy weak solution of ([I]) for any A G (0, 1). It was shown that for each initial data 
uq )X there is a unique, global strong measure-valued solution. The weak entropy solutions 
of (HJ are naturally contained in this approach as measures 

Vt,x = &u(t,x)- ( 7 ) 

The set strong measure-valued solutions is the subclass of the entropy measure-valued 
solutions were introduced by Tartar[10j in his compensated compactness method. Such 
solutions were further studied by DiPerna[4], who showed that the weak entropy solutions 
are unique in the class of measure- valued solutions and by Schochet[9], who showed that 
entropy measure-valued solutions with a prescribed initial data are not unique. 

An equivalent formulation of a strong measure- valued solution was given in Panov[8j 
where it was shown that Y(t, x, v) is the distribution of a strong measure valued solution 
iff it verifies ([5]). Note that the formulation (J5]) in addition to (J4|) prescribes a non-trivial, 
non- linear constraint: 

J J Y%x,v) dxdv = const., a.e.t>0. (8) 

This result was later re-discovered by Brenier[2] in the following form. He proved that Y 
is the solution of (|5|) iff u(t, x, A) from © is an entropy weak solution of (P) for any A in 
the interior of the range of Y. In that paper, the variational formulation ([5]) is expressed 
as a differential inclusion on a Hilbert space of I? integrable in (x, v) functions: 

Y t e — f v ■ V X Y — dK(Y), (9) 

where dK{Y) is the subdifferential to the indicator function of a convex, closed cone K 
consisting of all non-decreasing in v functions, see section[3]for details. Operator appearing 
on the right of © is monotone, providing the uniqueness of solutions and if in addition 
it is maximal then the existence follows from the classical results, for example Brezis[3i. 

Additionally to the existence/uniqueness/stability of solution Brenier[2j proves the reg- 
ularity of solutions of ([9]): if the initial data Yq(x,v) is differentiable, V^Yo £ ^ , then 

8 t Y, V X Y G L°°((0,+oo);L^„). 

It can also be shown that d v Y £ L°°((0, +oo); L x v ), if in addition 8 v Yq G l? x v . 

The results [2] show a remarkable fact that all weak entropy solutions of ([T]) can 
be obtained through ([6]) (or ([7])) from a globally stable and regular (if d x Yq G L x ) kinetic 
densities Y(t,x,v) - solutions of ([9]) (or ([5])). 

In this paper we further investigate the properties of solution of ©. Such solutions can 
be considered as curves Y(t) with values in the admissible cone K for which the tangent 
d t Y(t) belongs to the tangent cone T K (Y(t)) = closure {h(Y - Y(t)) : h > 0, Y G K}. 
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Our main result shows that Y(t) is the solution of ([9]) iff dfY(t) minimizes the functional 
II V + fv • V x y(i) 1 1 v over all directions V G Tx(Y(t)). This functional can considered as 
an "instantaneous" interaction functional. Since set K restrict solutions only in v direction 
this functional is local in x, i.e., 

min \\V + f v -V x Y(t)\\ L 2 = || min ||V(-) + /„(•) • V x Y(t, x, \\ L 2 

T K (Y(t)) J »^x,v T K (Y(t,x,-)) " 1 

Zero minimal value is attained on solutions Y(t) that are simply transported, i.e. Y t + 
fv ■ V X Y = 0, while remaining in K. At the level of weak entropy solutions of (pQ) such 
Y(t) corresponds to classical solutions. For shock waves it is proportional to the shock 
strength, see section El 

After the prove of this result, which based on the fact that solutions of 9(7 £ — 
A(Y(t)), with maximal monotone operator A are slow solutions, i.e., the solutions for 
which ||„4(y(i))|| is minimal, we show that there are travelling wave solution to Q that 
correspond to the shock waves of (JTJ) . Such travelling waves move with the actual shock 
speed a = A// Am. The shock speed appears in solving a minimization problem 

min \\V + f v -V x Y(t)\\ L 2 . 

The shock profiles of this type are obtained by smoothing in x direction the kinetic 
density (|3|) of the shock wave u(t,x). This however is rather exceptional case. In the last 
part of this paper we give an example that shows that generically if u(t, x) is a weak 
entropy solution that contain interacting waves and if Y(t, x, v) - its kinetic density, then 
Y £ = Y(t,x,v) * uj £ (x), is not a solutions of ©. This happens because the constraint (JSj) 
is non-linear in Y and does not commute with averaging. 

2. General theory 

Let H be the space of 2L-periodic in x, functions Y(x,u) of (x,u) £ W 1 x [0, 1], with 
the norm 

\\Y\\ 2 = (Y,Y) = I / Y 2 (x,u)dudx, U=[-L,L] n . 
in Jo 

Let ifcHba set of Y's non-decreasing in u. K is a closed cone and for any Y E K, we 
denote 

T K (Y) = M- closure of {h(Y -Y) : h>0,YeK}, (10) 

the tangent cone to K at Y and the normal cone: 

dK(Y) = {Z £ I : (Z,Y — Y) < 0, VY € K}. (11) 

Also by N = {Z G H : (Z,Y) < 0, MY G K} we denote the polar cone to K. We consider 
the Cauchy problem 

d t Y + f u -V x Y e -8K(Y), Y(t = 0) = Y . (12) 

The flux function / G Lip([0, l]) n - Lipschitz continuous on [0, 1]. 

It was shown in Brenier[2j that for any Yq G K, there is a unique solutions of (1121) 
Y G C([0, +oo); H), for any t > 0, Y(t) G K, and if V X Y G H, then also 

dtY, V x y G L°°((0,+oo);EI). 

Our main result contained in the following theorem. 
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Theorem 1. Let Y G K and V X Y G M. For the solution Y G C([0, +00); H) of (fT2l) . 
a.e. t > 0, 

l|a * y(t)l1 = z e a^(t)) l|Z + /u ' v * r WII' (13) 



Wt) + /«-V ie y(t)||= min ||V + /»-V s r(t)||. (14) 

VeT K (y(t)) 



and 



Conversely, each of the conditions f)13|) . (|14p defines a unique solution of the problem (|12p , 

Proo/. Let a(«) = /»/|/»|, « G [0,1]. If > 0, we define B(Y) = \f v (v)\d a{v) Y, 
where d a Y is the derivative of Y in the direction a. The domain of this operator, D(B) 
consists of all functions F6i that are weakly differentiable in the direction a(v) for a.e. 
v G [0, 1] and such that d a Y G H. 



Lemma 1. Let 



c = essinf \f v (v)\ > 0. (15) 



Then, B is a maximal monotone operator. 

Proof. Monotonicity of B follows directly from the definition of B and periodicity of Y in 
x. To show maximality, let W G H, and A > consider a problem: 

Y + X\f v \d a Y = W. 

For a.e. v G [0,1], Wo(v,-) G L 2 (U) and and the equation can be integrated along the 
characteristic to obtain a periodic solutions Y(v, ■). The inclusion Y G D(B) follows from 
the a priori estimates 

ll^ll < c ||d a y|| < (v^)' 1 !!^!!. 

□ 

Lemma 2. Under the condition on f from the previous lemma, B + dK is a maximal 
monotone operator on H. 



Proof. Consider now a proper, l.s.c, convex function 

K(Y) = 



0, x G K, 
+00, x g" K . 

The subdifferential dK is a maximal monotone operator. The Yosida approximation of 
dK(Y), equals VK X (Y), where 



K X (Y)= inf f K (Y) + ^-\\Y-Yf) 



see Theorem 4, p. 162 of Aubin-Cellina[I]. It follows that K\(Y) = infy^ ±\\Y - Y\\ 2 , 
and 

VK X (Y) = jTT N (Y), 

where ttn(Y) is a projection of Y onto iV - the polar cone to K. Operator jT^n(-) is a 
Lipschitz continuous operator from H to H, with the Lipschitz constant j and is monotone. 
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From lemma [I] and lemma 2.4 of Brezis[3] it follows that B + jTTn{ ) is a maximal 
monotone operator. So, for any W G H and a > 0, there is a solution Y x of 

Y + a(B(Y) + VK x (Y)) = W. (16) 



Using monotonicity we get 
Also, since 



|y A n < www. 



we obtain 

and consequently, 



(VK\Y x ),d a Y x ) = (±n N (Y x ),d a Y x ) = 0, 



a(\f v \d a Y x ,d a Y x ) = - (W,d a Y x ) 



Wd a Y x \\ <C(a,co)\\W\\, (17) 
||VK A (y A )|| < C(a,eBssvv\f v (v)\,co)\\W\\. (18) 

V 

We want to pass to the limit A — > in the equation (1160 . We have shown that all terms 
in that equation are weakly compact in M. It remains to show that the sequence Y x is 
strongly compact and use the strong-weak closeness of the maximal monotone operator 
dK. Using the equation (JIB]) we compute 

l|y A - y»\\ < a(\(VK x {Y x ), vF(y)) + /i(vf A (y A ),vF(y)) 

- A||V^ A (Y A )|| 2 - HlVi^Y^H 2 . 

This estimate, due to f)18[) . implies that Y x is Cauchy and converges to some Y G H. 
Moreover (I+XdK)- 1 ^) = ir K (Y) = Y X -XVK X (Y X ) converges to Y. Since VK X (Y X ) G 
dK({I + A<MT) -1 (y A )), and dK is strongly-weakly closed, it follows that VK X (Y X ) -> 
9ivT(y), and V is the solution of 

y + a(S(y) + aftr(y)) = w, 

proving by this the maximality of B + <9-fT. □ 
Consider a Cauchy problem 

<9iY g - B(y) - a^r(y), y(t = o) = y . (19) 

Under the non-degeneracy condition f)15f) , /3 + <9iT is maximal monotone and the problem 



(|19p has a unique solution with the properties listed in the next theorem, see theorem 1, 
p. 142 of pQ. 

Theorem. Let Yq G D{B) fl -fT. There is a unique solution Y(t) of (fT9j) for t G [0, +00), 
with the following properties: Y(t) G D(B) Pi K, 

Y G C([0,T];H), VT > 0, d t Y, d a Y G L°°(0, +00; H). 

Moreover, 

(1) // V a y G M, i/ien for any t > 0, ||V«y(t)|| < ||^y ||. 

(2) c\y(-) is continuous from the right on [0, +00) and ||9 t y(f)|| < esssup n \f v (v )| H^YoH. 
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(3) For any t > 0, 

d t Y(t) = argmin ||y||. (20) 

Ye-\f v \d a Y(t)~dK(Y(t)) 

Let / G Lip([0, l]) n and f £ be sequence of Lipschitz continuous vector functions such 
that: (i) f £ -4 / in C([0, l]) n ; (ii) / Ei „ -> /„, a.e. u G (0,1); (hi) ll/e^lli^o,!)).* - uniformly 
bounded; (iv) for all e G (0, £o)> mi "[o,i] l/e,v| > 0- For each / £ and Y$ G D(B)(lK, there is a 
solution y e that solves (|19p with / e and verifies the conclusions of the cited above theorem. 
It follows from the same theorem and assumptions on f £ that norms ||y £ (i)||, ||£^y e (£)||, 
||V x y e (t)|| are uniformly bounded in (t,e) G [0, +oo) x (0,£o)- Moreover, by monotonicity 
we obtain: 

\\Y £l (t)-Y £2 (t)\\ <t(esssup|/ eii „-/ ea ,«|)||5 1B yo||, Vei,e 2 G (0,e ), (21) 

i.e. Y £ is compact in C([0, +oo); H). With this and using the fact that dK, as a maximal 
monotone operator is strongly- weakly closed, we obtain Y = lim Y £ - the solution of (|12p 
with d t Y,V x Y G L°°((0,+oo);H). 

Let us prove property (|13p . Consider operator A = B + dK. It is monotone and has a 
maximal extension A. Thus Y(t,x,u) is also a solution of dtY G — .4 00, and by ([50]) . 
||a t y(t)|| = min Ze ^ (yW) ||Z||, for a.e. i. Since 5 4 y G ^.(y(t)), and ([13]) follows. 

Now we can prove (|14p . For the tangent cone T^:(y), defined in the beginning of this 
section, we have that V G T K {Y) iff MZ G 5if (Y), {V, Z) < 0, i.e. T K (Y) is the polar cone 
to a convex closed cone dK(Y). Property ([13]) states that dtY equals (J — 7rgif)(— S(y)), 
where 7raK- is the projector onto dK(Y). This can be stated equivalently, that dtY is the 
projection of —f u ■ V X Y onto Tk(Y), or 

\\d t Y(t)+B(Y(t))\\= min ||F + B(y(t))||, 

V€T K {Y(t)) 

for a.e. t. □ 



3. Examples 

Consider a scalar conservation law ([1]) in one dimension with a convex flux function 
f(u). We prescribe the initial data 



uq(x) = 




x G [-L,0] U [L/2,L], 
x G (0,L/2), 



with u + > u . The weak entropy solution u(t, x) of ([T]) consists (for small times) of a 
shock wave propagating from x = with the speed 

o- = {f{u + )-f(u-))/{u + -u-) (22) 

and a rarefaction wave centred at x = L/2. The solution has this structure until the 
moment the shock wave collides with the r-wave. Let us choose a small e > and consider 
the kinetic formulation for this problem. We define 



(23) 
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Y(t,x,v) is the solution of the variational problem ([5]), but it is not differentiable in x, 
and so we can not test it in the interaction functional dHJ). We will approximate Y(t, x, v) 
by 

Y E (0, x, v) = Y(0, x, v) * uj £ (x), 

where u £ (x) is the standard (supported on [x—e, x+e], non-negative, unit mass) smoothing 
kernel. 

It can be verified that Y e (t,x, v) for all small t, for which the shock wave and the r- 
wave in u(t, x) are separated by the distance larder than 4e, is the solution of © (or 
((SJ)). Indeed, Y e verifies because a convolution in x with a non-negative kernel doesn't 
change the structure of that equation, and, moreover it can be checked by a computation 
that the conservation property flSJ) holds as well. The structure of Y £ is simple; it consists 
of a smoothed shock wave : for x G (at — 2e, at + 2s), it equals 

{1 v > u + , 

+ 2e to £ (y)dy v G [iT , u+] , 
v <u~, 

and the part that corresponds to the regularization of the rarefaction wave. Next we 
would like to find the value of the interaction functional. Let us fix time t = and for 
x G [— L,L], consider 

^^o^^ + ^^^^^^'^W- 1 ^- (24) 

In the next lemma we show that the minimal value of is zero when x is in the range of 
the r-wave and it is proportional to the shock strength \u + — u~\ for x in the range of the 
shock discontinuity. For the x's in the later case, dtY e + ad x Y £ = 0, where a from ([22)1 . 

Lemma 3. £e£ e < L/8. TTie minimizer Vq of (|24p equals 

( -ad x Y £ (0,x,v) x G (— 2e, 2e), u G [0, 1], 
F = 9tr e (0,x,w) = <^ -A(v)5 x .y £ (0,2;,r;) x G (-2e + L/2,L/2 + 2e),v G [0,1], (25) 
[ otherwise, 

with a from (122|> . T/ie minimal value is proportional to the strength of the shock wave 
|u + — u~|, /or x G (— 4e, 4e), and is /or oi/ier values of x. 

Proof. Assume that <r > 0. The other case it treated similarly. The approximate initial 
datum Y £ (0,x,v), for x G (— 2e, 2e), equals 

{1 v > u + , 

f*+ 2£ io £ (y)dy vE[u-,u+], 

v < u~, 

and for x G (-2s + L/2, L/2 + 2e), 

{1 v > u+, 

/i/2-2e w £ (2/) ^ v G [n" , n+] , 

w < u~. 
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d x Y £ in the corresponding ranges equals 

{0 v > u + , 

uj £ (-x) v G [u~,u + ], 

v < u~, 

for x G (— 2e, 2e), and 

{0 w > u + , 

-u £ (L/2 - x) v £ [u-,u+], 
w < u~, 

for x G (-2e + L/2,L/2 + 2e). Then, with 

F G T x (y £ (0,x,-)) = L 2 ((0,l))-closure of {h(Y(v) -Y £ (0,x,v)) : h > 0, d v Y > 0}, 

||^ + AW^.y £ (0,x,r;)|| 2 = [ U \V(v)\ 2 dv+ I' \V(v)-l\ 2 dv 



+ / l^ + AWa^O,^)! 2 ^. (26) 

Notice that, due to the fact that Y £ (0,x,v) takes only three values, all functions V(v), 
such that V = 0, for v G [0,u~), V(v) = 1, for v G (n + , 1], and V(v) is non-decreasing on 
belong to Tk(Y £ (0, x, •)). Thus, 

min lly + A^Y^O,^)!! 2 

min 1" \V(v) + f v (v)d x Y £ (0,x,v)\ 2 dv. 

V'(v)>0,ve[u-,u+]J u - 

For x G (-2e + L/2,L/2 + 2e) we can take V = oj £ (L/2 - x)f v (v), for -u G [n",n + ]. Such 
V gives zero value of the functional. For the shock discontinuity range x G (— 2e, 2e), 
because f v (v)d x Y £ (0,x,v) is non-decreasing in v, minimum will be achieved on constant 
functions V(y) = c : 



min / \V(v) + f v (v)d x Y £ (0,x,v)\ 2 dv 

V'(v)>0, ve[u~,u+] J u - 



c 



pa 
J u~ 



min/ |c + f v (v)d x Y £ (0,x,v)\ 2 dv. 



The later functional is minimized for c = — ad x Y £ (0,x,v). This establishes (|25|) . It is 
easily verified that with such minimizer the value of (I24p is proportional to \u + — u~\. 

□ 

Next we will show that the regularizations Y £ = Y * cr £ (x) of the kinetic density Y of a 
weak entropy solution u(t, x) are not, in general, solutions of the variational problems ([9]) 
(or ((SJ))- For that we consider a conservation law: 

ut + ((u-h 2 ) x = 0, (27) 
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with 2L periodic data 

r l, xe [-l,o], 

u (x) = < 0, x £ (0,1), 

I 1, xe [1,L). 

The corresponding entropy solution u(t, x) consists of a stationary shock wave at x = 0, 
and a r-wave centred at x = 1 that propagates to the left with speed 1. Moreover, the 
values of u(t,x) in the r-wave depend linearly on j. Let, as in the previous example, 
Y(t,x,v) be the kinetic function of u(x,t) and Y £ (t,x,v) = Y(t,x,v) *u £ (x), where for 
the definiteness we take u> e (x) to be smooth, non- negative function, compactly supported 
on [— 2e, 2e] and equal to 1 on [— e, e\. It was shown in the previous example that for 
small times Y £ is a solution of Q. Consider time t = 1 — the moment the r-wave reaches 
shock. Y £ (l, x,v) is a smoothing in x direction of the characteristic function of a triangle 
{(x,v) : x 6 (0,2), v £ (x, 1)}, and Y £ (0,x,v) is a smoothing in x direction of the 
characteristic function of a square {(x,v) : x G (0,1), v £ (0,1)}, For all small e, one 
directly verifies that 

/L rl i*L rl 

/ \Y £ (l,x,v)\ 2 dxdv > / / \Y £ (Q,x,v)\ 2 dxdv, 
-L JO J-L JO 

violating the conservation property ([8]). 
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